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ABSTRACT 

The  study  of  steady  fluid  flow  through  porous  media  leads  to  the 
consideration  of  elliptic  free  boundary  problems  in  the  unknown  function 
"piezometric  head".  An  integration  along  the  vertical  direction  transforms 
these  free  boundary  problems  into  variational  (or  quasi-var iational) 
inequalities,  which  are  easier  to  study  both  from  the  theoretical  and  from 
the  numerical  point  of  view. 
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SIGNIFICANCE  AND  EXPLANATION  f 

A typical  free  boundary  problem  involves  motion  of  a fluid  in  a reqion,  » 

the  boundary  of  the  fluid  beinq  unknown  and  having  to  be  determined  as  part  I 

of  the  solution.  Problems  of  this  type  involvinq  flow  of  fluids  in  porous 
media  have  important  applications  in  ground  water  management,  oil  reservoir 
technology,  and  seepage  of  water  through  dams.  Apart  from  certain  special 
cases  that  can  be  solved  analytically,  computational  methods  for  solving  such 
problems  have  usually  involved  a laborious  trial  and  error  procedure  - guess 
the  position  of  the  free  boundary,  and  improve  the  guess  by  iteration. 

In  1971,  Baiocchi  showed  that  many  of  these  problems  could  be  reformulated 
as  variational  inequalities.  The  use  of  variational  inequalities  (variational 
problems  where  the  solution  is  subject  to  inequality  constraints)  and  their 

( 

generalization,  quasi-var iational  inequalities,  has  led  to  profound  develop- 
ments, such  as  existence  and  uniqueness  proofs  for  several  classical  free 
boundary  problems,  as  well  as  for  many  new  problems,  and  also  very  effective 
numerical  methods.  The  advantage  of  the  method  from  a computational  point  of 
view  is  that  the  position  of  the  free  boundary  appears  automatically  in  the 

r 

calculation,  without  the  need  for  guessing  and  iterating.  These  developments 
are  surveyed  in  this  report. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


FREE  BOUNDARY  PROBLEMS  AND  VARIATIONAL  INEQUALITIES 
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C.  Baiocchi 

Part  1.  The  simplest  type  of  free  boundary  problem. 

Let  us  consider  the  following  problem: 

PROBLEM  A.  Given: 

(1.1)  D:  a bounded  smooth  domain  in  »n, 

(1.2)  f,  g,  i|i:  smooth  functions  on  D;  iji  < g on  3D. 

We  look  for  a couple  {0,  u}  such  that: 

(1.3)  0 is  an  open  smooth  subset  of  D, 

(1.4)  u is  a smooth  function  on  0, 

(1.5)  -Au  - f in  0, 

(1.6)  u = g on  3!1  n 3d, 

(1.7)  u “ i|i  and  Vu  = Vi|i  on  30  O D. 

problem) : on  a domain  0 whose 
with  too  many  conditions  (see  (1.7))  on 


the  outward  normal  derivative) . 
instead  of  just  one  function  g,  we 
prescribe  a couple  (g^g^  and  a partition  of  3D,  {3d,  SjD};  then  we  replace  (1.6)  by: 

3u  3g2 

(1.6')  u - gx  on  3jD  n 30;  — - on  SjD  O 0. 

Remark  1.3.  Under  suitable  assumptions  on  the  data,  the  maximum  principle  will  imply  that: 
(1.8)  u > i|i  in  0. 

If  we  denote  by  u the  function  defined  by: 

# “ 
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Problem  A is  the  sinplest  f.b.p.  (free  boundary 
boundary  is  partially  unknown  we  must  solve  a problem 
the  unknown  part  of  30. 

Remark  1.1.  Condition  (1.7)  is  equivalent  to: 

(1.7*)  u - i(i  and  on  30  Cl  d (-^  denotes 

Remark  1.2.  A very  similar  problem  is  the  one  where. 
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('u(x) 

for  x # 

^ , 

(1.9) 

u(xi 

' |J|(X) 

for  x f 

D \ ft  , 

we  can  reformulate 

problem 

A in  the  following  form: 

PROBLEM  B . Given 

D,  f,  q, 

,4*  as  in  ( 1 . 1)  , 

(1.2) , we 

look  for 

(1.10) 

u is  smooth  in 

D ; u >_  i|i 

in  D 

(1.11) 

u = q 

on  3D  , 

y Setting  ft  = { x 

c D|  u(x) 

> iMx) } 

(1.12) 

; have  -Au  = f 

in  f. 

such  that 


As  a physical  situation  which  leads  to  mathematical  problems  of  this  tyne  let  us  consider 
an  elastic  membrane  above  a region  D C which  is  stretched  along  8D  at  the  height 
z = g(x,y) , subjected  to  a system  of  vertical  forces  F = f(x,y)  and  forced  to  stay  above  an  ob- 
stacle z = ii(x,y)t  let  C be  the  "contact  region"  and  ft  be  the  complement  of  C in  D,  say 
ft  = D \ C;  then,  assuming  that  the  membrane  is  homogeneous  and  with  unitary  elastic  coeffi- 
cient, the  shape  z = u(x,y)  of  the  membrane  at  the  eauilibrium  position  must  satisfy  problem 
A (or  problem  B) . 

Remark  1.4.  Problem  A,  as  well  as  problem  B,  does  not  contain  all  physical  conditions  which 
must  be  imposed;  we  will  come  back  to  that  point  in  a moment. 

Let  us  consider  the  analogous  one-dimensional  problem  with  the  simplest  assumptions  on 
the  data;  say  the  problem  of  an  elastic  string  of  length  fc,  fixed  at  its  endpoints,  stretched 
upon  an  obstacle  z = *Mx)  (the  string  is  supposed  weightless  and  no  forces  act  on  it) ; see 
Figure  1.1  for  which  D = ]0,£[  and  ft  = )0,8l  U 1y»£[» 
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Figure  1.1:  An  elastic  string  stretched  over  an  obstacle 


In  the  picture  is  drawn  the  true  physical  solution,  which  can  be  obtained  geometrical ly  bv 


drawing  the  straight  lines  from  (0,0)  and  (t,0)  tangent  to  s • + (x)i  however  it  is  easy  to 


see  that,  if  ij ) meets  horizontally  the  x-axis  at  the  mint  a 


tion  of  problem  B is  given  by  u^lx) 


the  set  SI  being  now  )0,a[  U )y,f| 


In  order  to  see  what  type  of  physical  condition  we  neqlv  «»-i  it  a *«e  » 


tially  three  ways,  which  lead  to  three  equivalent  iand  well 


1st  way:  Minimum  energy  principle.  Let  K be  the  »et  I • 


finite  elastic  energy"  in  the  Sobolev  space  H (r>) 


and  let  J(v)  be  the  energy  associated  with  v: 


then  the  equilibrium  shape  u must  satisfyi 


and  (1.15)  is  a well  posed  problem  (minimum  of  a quadratic  coercive  functional  on  a non 


we  come  back  to  the  general  two-dimensional  probli 


2nd  way : Virtual  work  principle.  Let  us  define,  for  any  u,v: 

<1.16)  a(u,v)  = //  Vu*Vv  dxdy  ; L(v)  = JJ  f v dxdy 

D D 

Remark  that  a(u,  u-v)  and  -L(u-v)  are  the  work  of  the  elastic  forces  (resp.  of  the  forces 
f)  in  order  to  pass  from  the  configuration  u to  the  configuration  v;  then  the  equilibrium 
shape  u must  satisfy: 

(1.17)  u c K;  a(u,  u-v)  <_  L(u-v)  V v f K 

(1.17)  is  a variational  inequality,  which  has  a unique  solution  ( ST AMPACCH I A 1 s theorem,  see 
[38];  in  our  present  case  the  result  is  quite  obvious  because  of  the  symmetrv  of  a,  namely 
a(u,v)  = a(v,u)  v u,v) ; we  can  remark  that  (1.17)  is  the  EULER  inequality  of  (1.15),  so  that 
problems  (1.15)  and  (1.17)  are  equivalent. 

3rd  way:  Balance  of  forces  and  reactions.  Let  us  denote  by  (j(x,y)  the  reaction  of  the  ob- 
stacle: we  must  have  u = -Au  - f,  so  that: 

(1.18)  -Au  >_  f in  D , 

(because  u must  be  upward  directed)  and 

(1.18)  support(-Au  - f)  C { (x,y)  e d|  u(x,y)  = i>(x,y)}  , 

(because  we  can  have  reaction  just  in  the  contact  region);  to  (1.18),  (1.19)  we  must  obviously 
add  u f K.  The  problem  in  this  formulation  is  often  written: 

(1.20)  u e H^(D)  ; u = g on  3D  ; 

(1.21)  u >_  \Ji  ; -Au  f ; (u-ij/)  (-Au  - f)  = 0 in  D ; 

but  this  is  a heuristic  formula,  the  product  (u-\|j)  (-Au  - f)  having  no  meaning  in  general;  on 

the  contrary  u e K and  (1.18),  (1.19)  is  a well  posed  problem  (it  is  very  easy  to  show  that 

it  is  equivalent  both  to  (1.15)  and  to  (1.17)). 

Remark  1,5.  From  (1.18),  (1.19)  we  get  opviously  (1.12);  the  reverse  implication  is  in  qeneral 
false,  as  we  showed  by  the  one-dimensional  counterexamt le;  in  particular  it  wan  the  strength- 
ened form  (1.18),  (1.19)  of  (1.12)  which  was  needed  (compare  with  Remark  1.4). 

In  order  to  see  that  the  (unique)  solution  of  the  variational  inequality  (1.17)  solves 
problem  B we  still  need  a "regularity"  result:  in  (1.10)  we  did  not  specify  how  smooth  u 


must  be,  but  problem  A suggests  that  along  3f!  n D u and  ty  and  their  gradients  must  agree; 
now,  from  u >_  \|>  and  the  definition  of  SI  (see  (1.12))  it  is  clear  that,  if  ip  is  a c* 
function,  this  agreement  will  hold  if  for  instance  we  have: 

(1.22)  The  solution  u of  (1.17)  is  C1 

Actually,  a well  known  result  of  LEWY  and  STAMPACCHIA  (see  [32])  implies  that,  under  reasonable 
hypotheses  of  smoothness  in  (1.1),  (1.2),  the  solution  u of  (1.17)  has  second  derivatives 
in  L^(D)  for  any  q < +«;  in  particular  (1.22)  holds  and  (1.21)  has  an  obvious  meaning. 

Before  discussing  a more  difficult  f.b.p.  let  us  sum  up  the  results  of  this  section:  The 
LEWY-STAMPACCHIA  regularity  result  is  the  key  for  interpreting  the  solution  of  a variational 
inequality  as  the  solution  of  a f.b.p.;  this  was  already  pointed  out  in  [32],  where  the  f.b.p. 
solved  by  a variational  inequality  was  a problem  of  lubrication  (instead  of  the  problem  for 
the  elastic  membrane.) 

Historical  and  bibliographical  note.  Variational  inequalities  like  (1.17),  when  we  assume 
that  K is  a closed  convex  non  empty  subset  of  a Hilbert  space  H,  a(u,v)  = (u,v)  (i.e., 

the  scalar  product  in  H)  and  L(v)  = (f,v)H  with  f e H fixed,  give  rise  to  the  problem  of 
the  projection  of  f on  H;  and  the  generalizations  for  L a linear  continuous  functional 


on  H,  and  a(u,v)  bilinear  continuous  coercive  and 


:ic  on  H * H are  quite  obvious. 


The  generalization  to  the  case  of  non  symmetric  a(u,v)  (which  does  not  correspond  to  a 
minimum  problem)  is  due  to  STAMPACCHIA  [38]  (see  also  LIONS -STAMPACCHIA  [33]  for  parabolic  in- 
equalities) . The  possiblity  of  interpreting  the  solution  of  a variational  inequality  as  the 
solution  of  a f.b.p.  was  not  obvious,  and  in  fact  it  is  strictly  connected  to  the  regularity 
results;  this  idea  was  firstly  developed  by  LEWY-STAMPACCHIA  [32]  , and  since  then  a lot  of 
papers  have  been  dedicated  to  the  argument;  see  e.g.,  [161,  [34]. 

Concerning  the  problem  of  regularity  for  free  boundaries,  (remark  that  we  did  not  dis- 
cuss the  smoothness  required  in  (1.3))  in  recent  years  great  progress  has  been  made;  see  e.g., 
[18],  [29],  [30]. 

For  further  details  on  variational  inequalities  one  can  see  the  books  [6]  (recently 


appeared)  and  [31]  (which  will  appear  soon) . 
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Part  2.  A f.b.p.  in  hydraulics. 

Let  us  consider  the  problem  of  a liquid  flowing  through  a porous  dam  under  the  following 
schematic  assumptions:  The  dam  is  homogeneous,  isotropic,  incompressible , with  vertical 
plane  parallel  walls,  on  a horizontal  impervious  base;  the  liquid  is  incompressible;  the  flow 
is  irrotational  and  steady  (in  particular  the  heights  of  the  reservoirs  in  the  picture  are 
constant) . 


Figure  2.1:  Porous  flow  through  a rectangular  dam. 

Figure  2.1  shows  the  cross  section  of  the  dam; 

D = 10, a[  x ]0,b[  is  the  dam, 

!5  - { (x,y)  £ d|  y < *(x)}  is  the  flow  region  . 

We  remark  that  ■rtx)  (and  in  particular  the  height  of  the  point  are  unknown!). 

Let  us  denote  by  p(x,y)  the  pressure  at  the  (x,y)  point  of  D (the  atmospheric 
pressure  being  measured  by  zero);  DARCY's  law  implies  that  (assuming  unitary  some  physical 
coefficients)  if  we  set: 

(2.1)  u(x,y)  = p(x,y)  y in  !2  , 

we  have 

5* 

(2.2)  velocity  = -Vu  in  fi  , 

(we  denote  by  (1  the  interior  of  the  wetted  region)  so  that  we  get  for  u the  following 
relations: 
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Au  = 0 in  n . 


(because  of  (2.2)  and  the  incompressibility) ; 


= 0 along  OA  and  along  {y  = v>(x)  } 


(because  they  are  streamlines) , 


u = y^  along  OB;  u = y2  along  AC  , 


(we  assume  hydrostatic  pressure  on  OB,  AC) ; 

(2.6)  u = y along  {y  = ^(x)  ) and  along  CC^  , 

(the  pressure  is  atmospheric) . 

It  is  obviously  a f.b.p.:  along  the  line  y = 'f(x)  (which  is  the  unknown  part  of  30) 
we  impose  both  Dirichlet  and  Neumann  conditions  (compare  with  (1.7*)  and  Remark  1.1;  see  also 

Remark  1.2  about  the  condition  on  OA) ; the  unique  (but  very  important)  difference  with  respect 

3u 

to  Problem  A is  in  the  fact  that  the  Neumann  condition  on  the  free  boundary  is  = 0 in- 
stead of  (remark  that  now  we  must  choose  i|>(x,y)  = y,  in  order  to  have  u = i|i  on 

30  n D) ' ' . 

Remark  2.1.  By  assuming  a little  regularity  on  u (e.g.,  u e CQ(0))  we  easily  get  from  the 
maximum  principle; 

(2.7)  u(x,y)  > y in  0 , 

(compare  with  (1.8));  on  the  other  hand  (2.7)  means  that  the  pressure  must  be  strictly  positive 
in  0,  which  is  physically  obvious. 

Let  us  reformulate  the  problem  in  terms  of  the  pressure  p (recall  that  p is  defined 
on  the  whole  of  D,  and  vanishes  outside  0);  we  will  assume  the  continuity  of  p (which  is 
physically  obvious)  and  the  fact  (also  obvious  from  the  physical  point  of  view)  that  the  flow 
has  a finite  kinetic  energy  (see  (2.9)  below) . 

Problem  C.  Given  a,  b,  y^,  y^  with; 


Remark  (with  notations  of  Remark  1.2)  that  we  could  choose  3^D  = OA;  32D  = 3D\0A; 
g^lx.y)  = maxly^.y)  on  OB;  g^x.y)  » b on  EF;  g^x.y)  - max(y2#y)  on  AF;  and  g2  = 0. 
We  have  i|>  g on  3^D  (instead  of  i(i  < g)  but  this  does  not  give  additional  problems. 


■ 

< 


(2.8)  a •"  0;  0 < y2  < /j  ' b , 

we  look  for  a function  p such  that,  settinq  D = 10, a[  x 10, b[,  and  denoting  the  spaces  of 
continuous  (continuously  differentiable)  functions  on  D vanishing  on  ID  by  Cq ( D)  (C^(D)), 

(2.9)  p < C (D)  O Hl(D)  , 

o 

(2.10)  p(0,y)  = (y.-y)+;  p(a,y)  = (y,-y)  + ; p(b,x)  = 0 , 


p(x,y)  ^ 0 in  D 


p(x,y)  > 0,  yt  e 10, y[  imply  p(x,yt)  > 0 


and  such  that,  setting: 


ft  = { (x,y)  c d|  p(x,y)  > 0}  , 


| 1/  V(P+y) *Vv  dxdy  = 0 , 

(2.14)  < a 

l V v e C1  ( 1 0 , a [ x It) 
o 

Remark  2.2.  (2.14)  Is  the  usual  weak  formulation  for: 

(2.15)  A(p+y)  = 0 in  ft;  =0  on  3ft  D (D  U OA)  . 

Remark  2.3.  The  meaning  of  (2.12)  is  that  the  set  ft  defined  through  (2.13)  (which  is  open 
thanks  to  (2.9))  must  have  the  form: 

(2.16)  ft  = { (x,y)  e d|  y < V»(x> } , 

for  a suitable  choice  of  remark  however  that  (2.12)  does  not  impose  any  regularity  pro- 
perty for  that  f ( •fi  must  a priori  be  just  a function  lower-semi-continuous  because  of 


We  have  already  pointed  out  that  the  double  condition  on  the  free  boundary  (p  = 0 and 

1 U* = 0)  is  different  from  the  one  we  encountered  in  problem  A and  B;  now  let  us  point 
3v 


This  hypothesis  is  a "regularity  assumption”  on  the  set  ft  defined  throuqh  (2.13)  above; 
see  also  Remark  2.3. 
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out  that,  if  we  ask  for  the  treatment  of  Problem  C under  the  variational  inequality  form,  the 


LEWY-STAMPACCHIA  Theorem  gives  a negative  answer:  if  the  function  p were  the  solution  of 

a variational  inequality,  we  would  have  = 0 (instead  of  ^ = 0)  along  {y  = ^(x)}. 

dV  d V 

An  other  way  of  seeing  this  is  the  fact  that  p cannot  have  (as  required  from  the  LEWY- 

STAMPACCHIA  Theorem)  second  derivatives  in  L^(D) ; in  fact,  from  (2.14),  choosing  v e C (D) , 

o 

we  get : 

C -Ap  = D x0*  in  the  distribution  sense  on  D , 

(2.17)  < 

• where  is  the  characteristic  function  of  S in  D 

00 

Let  us  give  the  (very  simple)  proof  of  (2.17).  Equation  (2.14)  with  v € Co(D)  gives  (we 
denote  by  < , > the  pairing  between  distributions  and  test  functions  on  D;  recall  that 
p 5 0 outside  0) : 

0 = //  V(p+y)*Vv  dxdy  = //  7p*Vv  dxdy  + //  Vy-7v  dxdy  = 

n n si 

= //  7p*Vv  dxdy  + |J  v dxdy  = < -Ap,v)  + <Xjj,  v > = 

D D y y 

= < -Ap,  v > - <DyXn,  v>  » < -Ap  - DyXn,  v>  . 

Formula  (2.17)  suggests  taking  an  integral  of  p with  respect  to  y as  a new  unknown 

function,  in  order  to  destroy  the  "bad"  operator  Dy  acting  on  X^J  so  let  us  define: 

b _ 

(2.18)  w(x,y)  = / p(x,t) dt  , (x,y)  e D . 

y 


Easy  computations  give: 


(2.19)  Aw  » Xn  in  D , 

b b 

w(x,b)  - 0;  w(0,y)  - / (y^— t)  dt:  w(a,y)  = J (y2~t)  dt  , 

y y 

w « -1  on  AB  w = Aw  - w = 0 on  AB  , 
yy  XX  yy 

so  that,  if  we  denote  by  g(x,y)  any  (smooth)  function  such  that  g(x,b)  = 0, 
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3(0, y)  = f (y  1 - t)+dt,  q(a,y)  = / (y 2 - t)+dt,  g(x,0)  linear  (so  that  g(x,0) 

y v 


v->  y,  -y. 


(a-x) ) we  get 


w = g on  3D 


Now,  recalling  (2.12),  (2.13)  and  p ^ 0,  we  get: 


w > 0 on  D , 


Q = { (x,v)  e d|  w(x,y)  > 0}  , 


so  that  we  can  replace  (2.19)  by: 


-Aw  > -1;  w(-Aw  -1)  = 0 on  D 


Now  (2.21),  (2.22)  look  like  (1.21)  with  . 5 0,  f = -1;  because  of  the  regularity  of  g 
(whose  first  derivatives  are  Lipshitz  continuous  on  each  edge  of  D)  we  can  conclude: 
Theorem  2.1.  There  exists  a uni-rue  w r1  ( D)  , with  second  derivatives  in  L^(D)  vq  < +“, 


such  that  (2.20),  (2.21),  '2. 


This  theorem  implies  that  Problem  C has  at  most  one  solution;  if  such  a solution  exists, 
it  can  be  evaluated  starting  from  v and  setting  SI  defined  through  (2.22)  and  p defined 
through  (compare  with  (2.18)); 


(2.24)  p(x,y)  = w)  (x,y)  Y(x,y)  c D ; 

the  main  difficulty  is  the  fact  that  it  is  not  obvious  that  0 defined  through  (2.22)  has 


the  form  (2.16).  Actually  we  can  prove  that  S)  defined  through  (2.22)  has  the  form  (2.16) 

and  moreover  the  corresponding  V is  very  smooth;  precisely: 

( •f  £ C ([0,aJ);  ■/>  is  strictly  decreasing;  f is 

(2.25)  ( 

^analytic  in  ]0,a(;  ^(0)  = y^;  *’(0)  = 0;  ^(a)  > y^,  f (a)  = -“ 
so  that  we  can  conclude: 

Theorem  2.2.  Problem  c has  a unique  solution;  the  upper  boundary  of  S)  has  the  further  re- 


gularity properties  listed  in  (2  25). 
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Remark  2.4.  If  we  ask  for  further  regularity  properties  of  p (which  a priori  satisfies 


just  (2.9))  we  find: 


(2.26) 


p c L "(D)  P„  € L (D)  , 

x y 


and  such  regularity  is  '' optimal”  because  actually  we  have: 
(2.27) 


d / L (D)  ; p i C (D)  . 

x y o 


Now  let  us  point  out  a typical  computational  difficulty  which  always  appears  when  a 
f.b.p.  is  solved  by  means  of  variational  inequalities.  Many  papers  and  books  deal  with  the 
numerical  solution  of  variational  inequalities  (e.g.,  see  (27));  so  that  we  can  assume  that 


we  have  constructed  a family  {w,  (x,y)  }.  „ of  functions  which,  as  h -*■  0,  converges  in  some 

n n>u 


sense,  to  the  solution  w(x,v)  of  (2.20),  (2.21),  (2.23);  setting  e.g.,  ph<x,y)  = -DyWh<x,y) 


(compare  with  (2.24))  we  will  also  have  ph  -*•  p in  some  sense;  but  how  can  we  approximate 
SI  (which  is  the  true  unknown  of  our  problem)?  The  naive  idea  of  defining 


S2.  = {(x,y)  £ d|  w.  (x,v)  > 0)  does  not  work,  no  matter  in  what  topology  w.  •*  w:  an  obvious 
n n n 


counter-example  is  given  by  choosing  wh<x,y)  = w(x,y)  + h (so  that  wh  -*■  w in  the  C 


topology!)  so  that  ft.  E D.  However  we  can  bypass  this  difficulty  by  means  of  the  following 
h 


theorem. 


Theorem  2.3.  Let  w.  w uniformly;  and  let  £,  be  such  that 
h h 


(2.28) 


V 0 • 


'l  (d) 


-*•  0 . 


Setting  £1  = {(x,y)  e dI  w (x,y)  > e ) we  have  "ft  -*•  ft  from  the  interior",  say  ft  = lim  ft 


h*0 


(lim  in  the  set  theoretic  sense)  and  ft.  C fi  for  h sufficiently  small. 

h 


Historical  and  bibliographical  note.  Problems  like  Problem  B (perhaps  in  a less  precise  form) 
are  very  ancient;  aoart  from  some  analog  solutions  (HELE-SHAW  model)  they  were  solved  numeri- 
cally, before  the  advent  of  computers,  by  means  of  the  so  called  "inverse  method"  or  "hodo- 
graphic  method";  see  e.g.,  the  book  (36).  For  general  problems  in  hydraulics  see  e.g.,  the 
books  [13],  (28),  (36).  After  the  advent  of  computers  a more  efficient  numerical  method  has 
been  proposed;  the  idea  is  the  following  one:  we  fix  an  initial  guess,  y « ^qIx) , for  the 
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free  boundary;  in  the  corresponding  set  we  solve  the  problem  (2.14)  and  the  two  first 

relations  of  (2.10);  the  solution  oQ(x,y)  will  not  satisfy  the  condition  nQ(x,  ^(x))  = 0 

' 

(otherwise  we  have  already  solved  our  problem!)  but  we  can  use  the  values  of  D^(x,  ^(x)) 

in  order  to  modify  (x)  and  get  a new  guess  <fi  (x)  (e.g.,  we  could  define 

.01  t 

^(x)  = V>0(x)  + Pq(x>  ^(x)));  a fixed  point  for  the  transformation  •+  will  give  the 

solution  of  our  problem.  This  is  just  a heuristic  approach  (e.g.,  with  the  formula 
^(x)  = '(’q(x)  + Pq(x»  ( x) ) we  will  never  modify  the  starting  value  for  \f(a)  l)  but  numeri- 
cally it  works  (for  references  to  this  method,  and  others,  see  [21)). 

The  transform  (2.18)  was  suggested  in  [2);  beside  the  existence  and  uniqueness  result  for 
the  problem  it  suggested  a new  numerical  treatment  (see  [20],  [8))  which  is  theoretically 
correct  and  which  competes  very  well  with  the  previous  heuristic  methods  both  in  simplicity 
of  programming  and  in  speed  of  execution. 

The  problem  of  the  convergence  ■*  Si  was  firstly  solved,  in  a particular  case,  in  ? 

[Ill;  the  general  result  of  Th.  2.3  has  been  exploited  in  [12],  in  connection  with  a parabolic 

f.b.p.  (diffusion-absorption  of  oxygen  in  tissues).  Remark  that  (2.28)  requires  the  knowledge 

00 

of  (or  an  estimate  for)  the  L norm  for  the  error;  such  an  estimate  (the  best  one: 

II w - w.  ||  = a(h^  E)  V£  > 0,  h mesh-size,  for  piece-wise  linear  elements)  has  been  obtained 

h °° 

L o 1 

independently  in  [3]  and  [35].  Some  improvements  of  Theorem  2.3  (say  L or  H norms  in- 
stead of  L norm)  can  be  given  (H.  Brezis,  personal  talk,  unpublished).  Concerning  (2.25),  ^ 

(2.26),  (2.27)  see  [18],  [29],  [30].  For  other  related  results  see  part  3 as  well  as  [6]  and 
its  bibliography. 


Part  3.  Recent  progresses  and  Bibliography. 


Because  of  space  restriction  I will  give  a few  details  just  for  some  other  problems  of 
elliptic  type  related  to  hydraulics;  however  let  me  point  out  that  transforms  like  (2.18)  have 
been  adapted  to  many  different  problems;  e.g.,  for  fluids  dynamics  problems  (flow  past  a sym- 
metric aerofoil;  see  [15],  [16]),  for  heat-conduction  problems  (like  the  classical  Stefan  pro- 
blems; see  [22]),  for  non-steadv  hydraulics  problems,  which  differ  from  both 
elliptic  problems  and  parabolic  ones  (in  fact  the  problem  in  the  interior  involves  only 

i 

space-derivatives,  the  evolution  condition  being  just  on  the  free  boundary)  see  [41], 
and,  for  more  details  and  references,  see  the  talk  of  Professor  FRIEDMAN  [24]). 

For  a general  overview  on  the  class  of  f.b.p.  which,  by  means  of  a suitable  change  of 
unknown  functions  (like  (2.18)),  can  be  translated  into  variational  inequalities  see  [5]  and 
the  Appendix  4,  Vol.  II  of  [6]. 

Now  let  me  give  some  more  details  concerning  steady  problems  in  hydraulics:  i 

A)  Variable  permeability  can  be  taken  into  account  with  some  restrictions;  e.g.,  horizontal 
layers  as  well  as  vertical  layers  can  be  studied  by  this  method;  see  [8] , [9] , [14] . 

B)  Phenomena  like  evaporation,  or  partially  pervious  bases  are  studied  in  [23],  [37]. 

C)  Capillarity  offers  some  theoretical  difficulties;  the  corresponding  numerical  treatment 
seems  to  suggest  that  the  theory  must  be  correct.  See  [19]. 

D)  More  general  geometry  can  be  taken  into  account;  when  the  right  wall  is  vertical  we  get 
(instead  of  one  variational  inequality)  a family  of  variational  inequalities  depending  on 

a real  parameter  (the  discharge  of  the  dam) ; the  value  for  the  parameter  is  found  by  impos- 
ing a "regularity  condition”  (see  [8],  [9]).  (then  also  the  right  wall  is  sloping  the  pro-  . - 

blem  becomes  a quasi-variational  inequality,  i.e.,  a family  of  variational  inequalities 
depending  on  a functional  parameter  (see  [4]). 

E)  Three  dimensional  problems,  in  the  case  of  vertical  walls,  do  not  offer  new  difficulties; 
see  [39];  for  non  vertical  walls  the  corresDonding  ouasi-variational  inequality  has  been 
studied  in  [26]. 

F)  A different  way  of  treating  the  problem  in  n-dimensions  and  with  very  general  geometry  has 
been  proposed  in  [1])  the  method  is  based  on  the  search  for  the  minimum  of  kinetic  energy 
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in  a suitable  class  of  supersolutions;  the  fact  that  the  point  of  minimum  satisfies  the 
problem  is  proved  by  a "balayage"  which  uses  locally  the  transform  (2.18). 

G)  Another  approach,  which  leads  to  a nonlinear  equation  (instead  of  a variational  inequality) 
has  been  proposed  in  [42] . 

H)  Problems  of  fluid  flow  through  a channel  are  studied  in  140) . 

I)  Problems  for  coastal  aquifers  can  also  be  studied;  in  [8) , [9)  was  treated  the  problem  with 
vertical  cliffs;  [7),  (10)  treat  the  problem  with  general  geometry  and  in  the  presence  of 
wells  sunk  into  the  aquifer  and  pumping  water;  conditions  which  ensure  that  just  fresh  water 
is  pumped,  and  not  the  salt  water  which  infiltrates  under  the  aquifer,  are  also  discussed. 


I 

i 
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Free  boundary  problems;  Variational  inequalities;  Porous  flow. 


20.  AIS^ACT  (Continue  on  rovoroo  oldo  II  noeoooory  end  Idontily  by  block  number) 

The  study  of  steady  fluid  flow  through  porous  media  leads  to  the  considera- 
tion of  elliptic  free  boundary  problems  in  the  unknown  function  *lpiezometric 
head^T  An  integration  along  the  vertical  direction  transforms  these  free 
boundary  problems  into  variational  (or  quasi-variational)  inequalities,  which 
are  easier  to  study  both  from  the  theoretical  and  from  the  numerical  point  of 
view. 
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